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In this paper, we discuss the angular momentum distribution in the ground states 
of many-body systems interacting via a two-body random ensemble. Beginning with a 
few simple examples, a simple approach to predict P(/)'s, angular momenta / ground 
state (g.s.) probabilities, of a few solvable cases, such as fermions in a small single- j 
shell and d boson systems, is given. This method is generalized to predict P(/)'s of more 
complicated cases, such as even or odd number of fermions in a large single- j shell or a 
many-j shell, d-boson, sd-boson or sdg-boson systems, etc. By this method we are able 
to tell which interactions are essential to produce a sizable P(I) in a many-body system. 
The g.s. probability of maximum angular momentum I ma x is discussed. An argument on 
the microscopic foundation of our approach, and certain matrix elements which are useful 
to understand the observed regularities, are also given or addressed in detail. The low 
seniority chain of g.s. by using the same set of two-body interactions is confirmed but it 
is noted that contribution to the total g.s. probability beyond this chain may be more 
important for even fermions in a single- j shell. Preliminary results by taking a displaced 
two-body random ensemble are presented for the / g.s. probabilities. 
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1 INTRODUCTION 



Many-body systems interacting via a scalar two-body random matrix elements are 
expected to have eigenstates which are (nearly) random superpositions of Slater de- 
terminants. Among all the quantum numbers, the particle number n and the total 
angular momentum / of the states are the only good quantum numbers. Regulari- 
ties exhibited by finite many-body systems interacting via random matrix elements 
provide an excellent tool to study general features, which are independent of in- 
teractions, of many-body systems 0. Therefore, robust regularities of many-body 
systems, if there are any, are both interesting and important. 

Recently, Johnson, Bertsch, and Dean discovered [|| that the dominance of + 
ground state (0 g.s.) of even fermion systems can be obtained by using a two-body 
random ensemble (TBRE). Further studies showed that quasi-ordered spectra can 
be obtained by using a two-body random ensemble || f|, |5|. The g.s. dominance 
was soon confirmed in sd-boson systems ]5j, ||] . Therefore, this g.s. dominance in 
even nucleon systems and boson systems is robust and insensitive to the detailed 
statistical properties of the random ensemble, suggesting that the features of pair- 
ing arise from a very large ensemble of two-body matrix elements and might be 
independent of the specific character of the force. An understanding of this g.s. 
dominance is important, because this observation seems to be contrary to the tra- 
ditional assumption which was taken previously. For example, in nuclear physics 
the g.s. dominance in even-even nuclei is explained as a reflection of a strong 
pairing associated with a strong short-range attraction between identical nucleons. 
Very recently, interesting studies are performed to check whether the spectroscopy 
with random and/or displaced random ensembles can simulate the realistic systems 

There have been a few efforts to understand this observation. In Ref . || , it was 
indicated that there is a correspondence between a large distribution width of + 
states and the g.s. dominance. In Ref. ||, it was suggested that for a system 
of interacting bosons the probability that the ground state has a certain value of 
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the angular momentum is not really fixed by the full distribution of eigenvalues, 
but rather by that of the lowest one. In |10| , Mulhall et al. discussed the g.s. 
dominance fermions in a single- j shell by using geometric chaoticity and uniformly 
distributed random interactions. Kusnezov discussed the sp-boson systems by using 



random polynomials [lT|. In Ref. [12, 13], the authors discussed the g.s. domi- 
nance in sp- and s<i-boson systems in terms of the mean field approach. In Ref. [|14| 
Kaplan et al. studied the correlation between eigenvalues and spins of corresponding 
states in a few simple cases. 

For fermions in a small single- j shell, it was recently shown [15| that the width of 
the energy distribution for each angular momentum I is not the key to understand 
the + g.s. dominance, instead, the coefficients ajppi (n-body matrix elements of 
two-body interaction A J ^ ■ A J ) with (3 = (3' was suggested to provide a reasonable 
explanation of the distribution of angular momentum I g.s., where otjppi = (nf3I\ 
A Jjl ■ A J \nf3'I) with n being the particle number, (3 additional quantum number 
necessary to label the state, and I being the angular momentum of the state. The 
A Jj; ■ A J will be defined later. It was assumed [|l^, |16| that the off-diagonal matrix 
elements ajppi {(3 ^ f3') are neglected as an approximation, which is based on an 
observation |T5| that they are very small compared with diagonal matrix elements 



ajpp for / > states of even fermions in a small single- j shell, e.g., j = |, -ji if one 
uses a seniority conserved basis [IT]. This method has a disadvantage that it is not 
applicable to more complicated cases, e.g., fermions in a large single- j shell, where 
the off-diagonal matrix elements become also important. 

These studies are interesting and important, and have potentially impacted our 
understanding on the origin of one of the most characteristic features of nuclear 
spectra. All these approaches, however, address only simple or very specific (sp and 
sd bosons, or fermions in a small single-j shell) cases. None of these arguments 
could explain the regularities of angular momenta / g.s. probabilities such as those 



observed in JL5) . It is therefore very desirable to construct a universal approach to 
understand the g.s. dominance of even nucleon systems and sd-boson systems, 
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and meanwhile, to understand all the / g.s. probabilities (denoted as P(I)) of very 
different systems including both even and odd number of fermions in a single-j or 
a many-j shell. In this paper, we shall present an approach to predict the P(/)'s of 
all types of systems in a simple but universal procedure. 

In this paper, we use Gj's be a set of Gaussian-type random numbers with a 
width being 1 and an average being 0: 

p(Gj) = -^exp(-G 2 J /2), J = 0, 2, ■ • • , 2j - 1, (1) 

where the Gj's define the two-body matrix elements of fermions in a single-j shell 
as follows: 

H = J2 GjA J] ■ A J = J2 V2J+lGj (A J t x A J )° , 
j j 

A J t = -^(a}xo}) J , A^-^fexa/, Gj = (j 2 J\V\fj). 

For fermions in a many-j shell we use Gj^j^jzj^s in stead of G j's. For pure d 
boson systems, there are only 3 independent two-body matrix elements, parameter- 
ized by Co, C2 and C4 |18] which will be defined later. For s<i-boson systems we shall 
also define the two-body hamiltonian separately. All the results, except a few cases 
which we shall note clearly, are obtained by taking the two-body matrix elements 
to be a two-body random ensemble (TBRE) described by Eq. ([[]). 

This paper is organized as follows: In Sec. 2, We begin with a few simple but 

non-trivial cases a j — I shell with n = 3 and 4, and d-boson systems with the 

number of bosons running from 3 to 48. In these cases, eigenvalues are analytically 
given in terms of linear combinations of two-body matrix elements. These examples 
are very helpful to understand our simple approach proposed in this paper, and also 
interesting because they are analytically given. In Sections 3 and 4, we generalize 
our method and apply it to both even and odd numbers of fermions in a large 
single-j or a two-j shell, sd- and sdg-boson systems, where this approach continues 
to work reasonably well. For fermions in a single-j shell it is found in this paper 
that the G2 = — 1 always produces I = n ground state for even n fermions and 
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I — j — (n — l)/2 ground state for odd n. In Sec. 5 we discuss the foundation of 
our simple approach, and suggest that the essential part of the I g.s. probability 
is related to interaction Gj's which produce angular momentum I g.s. if Gj = —1 
and others are zero. This disproves a popular idea that the g.s. dominance of even 
fermions and bosons comes intrinsically from the two-body nature of the interactions 
and might be independent of the two-body interactions that used. In Sec. 6, we 
address the I ma x g-s. probability. For fermions in a single- j shell. We present the 
eigenvalues of the I max states with n =3, 4, 5, 6. This picture can be generalized to 
explain a few features of the I max g.s. probabilities of fermions in a two-j shell. In 
Sec. 7, We shall check the previous statements of "pairing" and seniority "chain" 
related to the g.s. dominance for fermions in a single- j shell, where seniority is well 
defined. We confirm the finding of a low seniority chain suggested in but note that 
the contribution to the g.s. probability beyond this chain may be more important. 
A summary of this work will be given in Sec. 8. In appendix A we present our 
preliminary results by using a displaced TBRE hamiltonian. In appendix B we give 
a simple algorism to calculate diagonal matrix elements of fermions in a single- j 
shell. In Appendix C, a few counter examples of the g.s. dominance, and counter 
examples of I = j for odd fermins systems are also given. 
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2 A FEW SIMPLE SYSTEMS 

Generally speaking, eigenvalues of a many-body system are not linear in terms of 
two-body matrix elements. In some simple cases such as fermions in a single- j shell 
with j < 7/2 and pure c/-boson systems, however, all eigenvalues are given in terms of 
linear combinations of two-body matrix elements. For these cases an understanding 



of the I g.s. probabilities using a concept of shift described in Ref. [13 is applicable: 
A state which has the largest and/or the smallest coefficients for a given Gj in the 
eigenvalues is favored to be the ground state. These examples will provide a useful 
clue to obtain a universal approach. 

2.1 j = | and | shells 

Fermions in a single- j (j = | or ~) shell are simple but non-trivial cases. An 
analytical relation between eigen-energies Ej and two-body matrix elements Gj are 
available. 

The eigen-energies of states with 3 fermions in a j — | shell are given by 

f 15 r + 6 g 

E5/2 = —Gq + -G 2 + 2^ 4 ' 

E 9/2 = JJ G2 + li^ 4 ' ^ 

The eigen-energies Ej of states with 3 fermions in the 3 = \ shell are as follows: 

^3/2 — 14^*2 + ff G4 + OGq, 

rp XX r 1 1 2 r 1 1 65/^ 

-£^5/2 — "6" Lr 2 + 77*^4 + 66^65 

-E7/2 — f Go + Y^^2 + + Y§G6, 

771 _ 13^ , 150 ^-y , mr< 

-Z9/2 — 42^2 ' 77" ( - J 4 + g^e; 

5/^ 1 13 n 1 52, 



-C/ll/2 — gLf2 + 22^4 ' 33^6' 

-E'l5/2 — #G*2 + M^ 4 M^ 6 ' (3) 
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For four fermions in a j — I shell, each eigenstate is labeled using seniority 
number v and angular momentum /. The eigen-energies Ej^ of states with 4 
fermions follows: 
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In Eqs. bold (italic) font is used for coefficients which are the largest (smallest) 

among all the I(v) states for a given J. 
We rewrite Eqs. (@-§) as follows: 

Eiiv) = E a Hv)Gj = £ a{Gj, (5) 



where k = I (3 (3 (—I(v) in this subsection), and a{ satisfies a sum-rule |L7 



XX = in(n-l). (6) 
j z 

A method to calculate the above a( for 4 fermions in a single- j shell is given in 
terms of 9-j coefficients in appendix A, where explicit expressions are available for 
v — states. The ai max of fermions in a single- j shell with n = 3 to 6 is given in 
sec. 2.4. 

By using a TBRE hamiltonian described by Eq. (|l|) and the eigen-energies given 
by Eqs. (0-§), it is easy to obtain the probability for each / ground state (J g.s.). 
On the other hand, one can predict the / g.s. probability without running a TBRE 
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hamiltonian. For example, the exact g.s. probability of 4 fermions in the j ' — I 
shell is determined by the following integral: 

dG I dG 2 j dG 4 f dG e f dE f dE 2{2) ■■■ f dE 8 

* (^o(o) - E a J m Gjj ■ ■ ■ 6 \E 8{4) - £ a J m Gj^j p(G )p(G 2 )p(G 4 )p(G 6 ). (7) 

The P(/)'s of fermions in a single- j shell with j = |, or | are given in Tables I- 
III. The row "TBRE" corresponds to results obtained by using Eqs. (0-f|) and 1000 
sets of a TBRE hamiltonian. The row "predl." corresponds to the probabilities 
calculated by integrals such as Eq. (^) for the g.s. probability of 4 fermions in a 
3 = \ shell. Using Eqs. the distribution width, gi( v p), of each state, is equal to 



J2j ( a /(t,)) • These widths are listed in the last row of Tables I-III. 
It is noticed that an understanding in terms of shifts works well: a state with 
one or more largest (or the smallest) oij, v \ has a very large probability to be the 
ground state or the highest state. The P(/)'s of states without the largest and/or 
the otj^ for a given J are very small. A schematic argument of such an observation 



was given in detail in [|15|, [L6| 



For 3 fermions in a j — | shell, all states have the largest otj. Therefore, the 
above argument by using shifts predicts that all P(J)'s are large, because all the 
3 states have the largest aj for a given J . In Table I, all states have indeed large 
probabilities, which are obtained by running 1000 sets of a TBRE hamiltonian and 
shown in the row "TBRE" , be the ground. The state with / = | have only one aj 
which is the largest among different / states, while the other two states with I — |, 
| have both the largest and the smallest aj. Therefore, P(I = |) is smaller than 
P(I = |) and P(J = |). 

For 3 fermions in a j — | shell, states with / = |, |, ^ have both the largest 
and smallest aj, and the state with I = \ has the largest aj (J = 0). Therefore, 
the above argument predicts that the probabilities of / = |, |, | and y g.s. are 
large. From Table II it is noticed that all these states have large probabilities to be 
the ground states, but, the state 1 = 1 state does not have a coefficient which is the 



smallest among a/'s with a given J but different /, indicating that the I — ~ g.s. 
probability is a bit smaller than those of I = |, |, y. The / = | and y states do 
not have coefficients a/ which are the largest or the smallest among all the states. 
The predicted probabilities of these two I g.s. are small. Note that it can be shown 
that the / = y cannot be the ground state. All these features are confirmed by the 
P(/)'s, which are calculated by running 1000 sets of a TBRE hamiltonian, in the 
column "TBRE". 

The situation is very similar in the case of 4 fermions in j — | shell: in each 
of the four states (I(v)= + , 2 + (4), 4 + (4), 8 + ) which have large probabilities as the 
ground state (and the highest state), there are coefficients a/^ which are the largest 
and/or the smallest for different I{v). Other states which do not have the largest 
and/or the smallest aj coefficients have very small probabilities to be the ground 
state or the highest state. Refer to the column "TBRE" in Table III. 

In Tables I-III, the P(/)'s obtained by running 1000 sets of a TBRE hamiltonian, 
and those obtained by using integrals such as Eq. (0) for the g.s. probability of 
4 fermions for a j — | shell, are well consistent with each other. Therefore, the / 
g.s. probabilities P{I) in the above cases are explained in terms of shifts, defined in 
Ref. 0. 

However, the P(/)'s in Eq. (|7]) is not yet within reach of a simple procedure by 
hand, and one has to evaluate this integral numerically. It is then very interesting 
and important to find a simple method to evaluate the / g.s. probabilities. 

Let M[ be the number of both the smallest and largest otj^ with a fixed J 
for a certain /, the / g.s. probability is approximately given by J\fj/{N m ), where 
N m = 2N — 1, N is the number of two-body matrix elements. The predicted I g.s. 
probabilities by this method are also given in the row "pred2." of Tables I-III. A 
reasonable agreement is obtained, though there are small differences, compared with 
those obtained by running a TBRE hamiltonian. Note that in the above examples, 
we use N m = 2N — 1 because all ^ 0)'s are (there is no smallest afc?). 
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2.2 (f-boson systems 

Similar to fermions in a small single- j shell (j — | or j — |), the relation between 
the two-body matrix elements and the eigenvalues of d-boson systems is also linear. 
The two-body hamiltonian of a o?-boson system is given by 

H d = ^MTld ((rf 1 " x d)) 1 x (d x d)y (8) 
From Eq. (2.79) of Ref.@, we have 

E = E + a'^n d {n d - 1) + [n d {n d + 3) - v(v + 3)] + j' [1(1 + 1) - 6n d ] , (9) 

where Eq contributes only to binding energies and not to excitation energies, n d is 
the number of d bosons. Eq. (|9]) can be rewritten below: 

E(v, n A , I) = E' (n d ) - (3'v(v + 3) + 7 7(I + 1). (10) 

We cite Eq. (2.82) of Ref. @: 

a' + 87' = c 4 , 
a' - 67' = c 2 , 
a' + 10/3' - 12 7 ' = c , 



which may be rewritten below: 

a' = i(4c 2 + 3c 4 ), 

= ^(7c -10c 2 + 3c 4 ), 

7=^(-c 2 + c 4 ). (11) 

Substituting these coefficients 7' into Eq. (|10|), and taking the two-body 
matrix elements Co,c 2 and c 4 to be the TBRE defined in Eq. (fj), one obtains 
that only 1=0, 2, and I max (= 2n) have sizable / g.s. probabilities (other / g.s. 
probabilities are zero), which are shown in Fig. 1. We notice following regularities 
of P(J)'s vs. n d . 
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1. The P(I m ax) is almost a constant (around 40%) for all n d (< 4); 

2. The -P(O) and P(2) are periodical, with a period 5(rirf)=6. 

3. All the P(Imax), P(0) and P(2) are near to 0, 20%, 40%, or 60%. The other 
P(/)'s are always zero. 

Below we explain these observations. From Eq. (p~0|) and Eq. (|TT|), we have 

c = 1, c 2 = c 4 = : E(v, n A , I) = E' (n d ) - j^ v ( v + 3 ); 

c 2 = 1, c = c 4 = : E(v, n A , I) = E' (n d ) + -v{v + 3) - -^1(1 + 1); 

c 4 = 1, c = c 2 = : E{v, n A , I) = E' (n d ) - ^v{v + 3) + ^1(1 + 1), (12) 

where E' (n d ) is a constant for all states. Based on Eq. (|T2"D, we obtain TABLE IV, 
which presents the angular momenta giving the largest (smallest) eigenvalues when 
q = — 1 (I =0, 2, 4) and other parameters are for d boson systems. In Table IV, k 
is a non- negative integer, and n d > 3. These angular momenta appear periodically, 
originating from the reduction rule of U(5)^SO(3). From TABLE IV, one notices 
again that a certain P(J) is large if one state with angular momentum I has the 
largest and/or the smallest a l I/3 (Eq. (H)) for a given / (I = 0,2,4). 

Note that when one searches for the smallest eigenvalue with Co = — 1 and 
c 2 = c 4 = in case A of Eq. (|T2|), one finds that many I states are degenerate at the 
lowest value. Therefore, again, we use N m = 3N — 1 in predicting the P(/)'s by the 
formula P(I) = Mi/N m . The results are well consistent, without any exceptions, 
with those obtained by running a TBRE hamiltonian. Take n = 4 case as an 
example, A/}=o = 3 and Mi=i max = 2. We predict that g.s. probability is 60% and 
Imax = 8 g.s. probability is 40% while all other / g.s. probabilities are zero. The 
g.s. and I max g-s. probabilities given by diagonalizing a TBRE hamiltonian are 
60.7% and 39.3%, respectively, and all other / g.s. probabilities are zero. Note that 
our predicted P{I)s of <i-boson systems are always consistent with those obtained 
by diagonalizing a TBRE hamiltonian in the examples that checked: n^=4 to 48. 
In another sentence, the distribution of the P(/)'s can be explained satisfactorily 
by shifts produced by the largest and/or smallest a l k . 
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3 FERMIONS IN A LARGE SINGLE- j SHELL 



In this subsection we generalize the method proposed above, and study the P(/)'s 
of fermions in a single- j shell with particle number n=4, 5, 6, 7, 8. The explanation 
of our approach will be addressed in Sec. 8. 

The procedures of our approach to a general case are as follows. First, set one of 
the two-body matrix elements to be —1 and all other interactions to be zero. Then 
one finds which angular momentum / gives the lowest eigenvalue among all the 
eigenvalues of the full shell model space. Suppose that the number of independent 
two-body matrix elements is N, then the above procedure is iterated N times. Each 
time only one of the G/s is set to be —1 while all the others are switched off. Next, 
among the N runs one counts how many times (denoted as Hi) of a certain angular 
momentum / gives the lowest eigenvalue among all the possible eigenvalues. Finally, 
the probability of / g.s. is given by Ni/N m x 100%. Below we use N m = N (unless 
pointed out explicitly) although the largest eigenvalues are equivalent to the lowest 
eigenvalues. The reason is that these largest eigenvalues are usually (exactly or 
nearly) zero for many / matrices, especially for a many-j shell or a large (j > 9/2) 
single- j shell. To have the "rule" as simple as possible, we shall use only the lowest 
eigenvalues with one of the Gj's being set to be —1 and others being switched off for 
fermion systems in a large single- j shell, a many-j shell, sd- and s<ig-boson systems. 

Fig. 2 presents a comparison between the predicted P(/)'s and those obtained 
by diagonalizing a TBRE hamiltonian of fermions in a j — | shell. We present 
two cases: n = 4 (even) and n = 5 (odd). The agreements are good. Note that 
such agreements do not deteriorate (or become even better) if one goes to cases of 
fermions in a larger single- j shell or a many-j shell where there are more two-body 
matrix elements. 

Tables V-IX give the angular momenta / which produce the lowest eigenvalues 
for different two-body matrix elements and particle numbers, with Gj being —1 and 
others being 0, and with n as much as possible. In Table VI, the number of Mq 
staggers with j at a period of Sj—3, which will certainly produce a staggering of the 
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g.s. probabilities with j for 4 fermion systems. Fig. 3 gives a comparison between 
the predicted P(0)'s (open squares) and those obtained by diagonalizing a TBRE 
hamiltonian (solid squares) for n = 4 and 6. It could be seen that a good agreement 
is obtained for fermions in both a small single- j shell and a large single- j shell. 
The predicted g.s. probabilities exhibit a similar staggering as those obtained by 
diagonalizing a TBRE hamiltonian. 

It is interesting to note that P(0)'s can also be fitted by empirical formulae. For 
example, P(0)'s can be predicted by 



[(2j + l)/6] + k 



1 if 2j = 3m 



for n = 4 : P(0) = LV , [\ 1 x 100%, k = { if 2j + 1 = 3m 

-1 if 2j - 1 = 3m 



2 



forn = 6 : P(0) = 12^1 x 100%, (13) 
3-2 

where the " [ ] " means to take the integer part. These empirical formulas are 
interesting because it presents a scenario without any calculations for very large- j 
cases where it would be too time-consuming to diagonalize a TBRE hamiltonian. 

A very interesting note of Tables V-IX is on the quadruple matrix elements G2 
term. It has been well known for a few decades, based on the seniority scheme, that 
the monopole pairing interaction always gives 1 = ground state for even fermion 
systems in a single- j shell and I = j ground state for odd number of fermions in 
a single-j shell when Go is set to be -1 and others 0. However, little was known 
about the G2 matrix elements in a single-j shell. The Tables V-IX show that the 
quadruple pairing interaction corresponding to G2 always gives I = n ground state 
for even fermion systems and I = j — (n — l)/2 ground state for an odd number of 
fermions when G 2 is set to be -1 and others 0. A study of this observation based on 



pair approximation is now in progress |T9] 
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4 sd and sdg BOSONS, AND FERMIONS IN A 
MANY-j SHELL 



Although sd- and s<ig-boson systems, and even and odd numbers of fermions in a 
many-j shell are very different systems from the cases discussed above, our method 
is applicable. All features are explained similarly 

The hamiltonian of a sd-boson system is as follows: 

H s d = H d + e M „^(sV)(ss) + e sddd (J^i^^) {dd^ + h. 

+e ss dd (^\J\( sfsj ) (^0° + h - c ^j + e *dsd ((s f rf f ) x (sd) 

where H d is a two-body hamiltonian defined in Eq. (|S]). 

TABLE X presents the angular momenta which give the lowest energies when 
one of the above parameters is set to be -1 and others 0. We predict, according to 
Table X, that only I = 0,2, In g.s. probabilities are sizable, which is consistent with 



~ \ 



the previous observation || [13j , that in s<i-boson systems interacting via a two-body 
random ensemble only / = 0, 2, and 2n (maximum) have large probabilities to be 
the ground state, the g.s. probabilities of other angular momenta are close to zero. 

Fig. 4 shows a comparison of the predicted -P(/)'s and those obtained by diag- 
onalizing a TBRE hamiltonian of s<i-boson systems, with boson numbers ranging 
from 6 to 16. It is seen that a good agreement is obtained. 

The case of fermions in a many-j shell is the most complicated. Fig. 5 present 
a detailed comparison of the predicted P(/)'s and those obtained by diagonalizing 
a TBRE hamiltonian, for fermions in a two-j (j = ~, |) shell with n=4 to 6. The 
predicted P(J)'s are reasonably consistent with those obtained by diagonalizing a 
TBRE hamiltonian. 

For fermions in a many-j shell, number of two-body matrix elements is usually 
large. In such cases, especially in odd-fermion systems, there are "quasi-degeneracy" 
problem in counting A/j: sometimes the lowest eigen-value is quite close to the 
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second lowest one when one uses Gj^j^j^ = — 1 and others 0. For such two-body 
matrix elements, one should actually introduce an additional "rule" in order to have 
a more reliable prediction. Namely it is not appropriate to count M% in a simple 
procedure. In order to avoid confusions, however, we did not modify the way in 
counting Mx of such cases throughout this paper. It is noted that the I = | in Fig. 
5b) and / = | in Fig. 5d) belong to the case with "quasi-degeneracy" . Improvement 
of agreement between the predicted P(/)'s and those obtained by diagonalizing a 
TBRE hamiltonian can be achieved by appropriately considering the above "quasi- 
degeneracy" . 

We have checked two-j shells such as (2ji, 2j 2 ) = (5, 7), (5, 9), (11, 3), (11,5), 
(11,9) and (13,9) with n = 4,5, 6, sci-boson systems with n up to 17, and scig-boson 
systems with n — 4,5, and 6, and all the agreements are reasonably good. 
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5 A SCHEMATIC INTERPRETATION OF OUR 
APPROACH 



It is interesting and very important to know why our simple approach by numerical 
experiments can successfully produce / g.s. probabilities which are in good agree- 
ment with those obtained by diagonalizing a TBRE hamiltonian. In this subsection, 
we shall provide a schematic explanation. A sound explanation may be much more 
sophisticated. 

As mentioned above, the relation between the eigenvalues and the two-body 
matrix elements is usually not linear. However, eigenvalues are always linear in 
terms of two-body matrix elements in a "local" space (explained below). Namely, 
within the local space we may find linear relations between the eigenvalues and 
the two-body matrix elements. Therefore, instead of studying the effects of all 
the two-body matrix elements simultaneously, we dismantle the problem into N 
parts. In each part we focus on only one term of two-body matrix elements. As 
a schematic interpretation of our method, we take n fermions in a single- j shell 
as an example (It is easily recoginzed that this explanation is applicable to all 
complicated cases as well). Let us take a certain Gj = — 1 and all Gj> — (J' ^ J), 
and diagonalize the two-body hamiltonian. Suppose that the eigenvalues are Ej^, 
and their corresponding wavefunctions are 

$(j'V/3J)= E U n -^)j 2 (K')\}j n I/3J)(^(j n -\K 7 )) x*(j*(K'))y. (14) 

KK'-y 

Now we introduce a small perturbation by adding {eGj>}. Gj = — 1 and {eGj/} 
define our (J/2)-th local space of two-body matrix elements. Then the new eigen- 
energies are approximated by 

K)' = ^i + e ^f^ e [(r-^fmiyripj^Gj, (is) 

Z KK'jJ' 

Namely, the Ef^s are linear in terms of {Gj/} in the local space. For two-body 
matrix elements which are close to the above local space, the angular momentum 
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/, which gives the lowest eigenvalue among all I n s with Gj = — 1 and others zero, 
continues to give the lowest eigenvalue. This means that the a very large part of full 
space of a TBRE hamiltonian can be covered by the N = j+\ local subspaces defined 
above, especially one uses a TBRE hamiltonian which produce a large probability 
for small |Gj|'s. This is the phenomenology of our approach to predict the / g.s. 
probabilities in this paper. 

A further rationale can be seen from the following analysis. To exemplify briefly, 
let us take 4 fermions in a single- j shell with j = y. In Fig. 6a) we set Gj max (J max = 
16) = —1, and set all the other parameters are taken to be the TBRE but with a 
factor e multiplied. We see that almost all cases of the g.s. belong to / = I max when 
e is small (say, 0.4). If one uses Gj max (J max = 16) = 1, then the P(I ma x) ~ 0, which 
means that the cases of the TBRE with G j max < produce almost all the I max g.s. 
in a single- j shell. In Fig. 6b) we present the results of the same system with Go 
being —1 and other Gj's being the TBRE multiplied by e. It is seen similarly that 
the g.s. is dominant for small e, and that if we switch off all the interactions which 
give the 1 = lowest eigenvalue with a certain Gj being —1 and others 0, then the 
g.s. probabilities will be very small (such as 10%) or be close to zero (~ 2%). 
Therefore, by this method one readily find which interactions, not only monopole 
pairing, are important to favor the g.s. dominance. Previously, Johnson et al. 
noticed that the g.s. dominance is even independent of monopole pairing 0, |^, M. 
It was not known, however, whether a certain two-body matrix element is essential 
or partly responsible, and how to find which interactions are essential, in producing 
the g.s. dominance for a given system. 

A shortcoming of the above explanation is as follows: in our simple approach we 
set each Gj = — 1 for each numerical experiment and find the angular momentum 
of the lowest state. In most cases we obtain degenerate lowest states if we set 
Gj = 1. Thus the local space of {Gj = 1 + eGj'(J' ^ J)} is not considered 
according to the above explanation. However, the good consistence of our predicted 
/ g.s. probabilities with those obtained by diagonalizing a TBRE hamiltonian seems 
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to indicate that the local space such as {Gj = 1 + eGj>(J' ^ J)} is considered 
via certain procedures, namely, the properties of local spaces defined by {Gj = 
— 1 + eGji(J' ^ J)} are, more or less, enough to represent the main features of the 
full space, suggesting that the total / g.s. probilities of Gj = — 1 local spaces might 
be symmetric as a whole to those of the Gj = 1 local spaces. 
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6 The I max g.s. PROBALITIES 

For fermions in a single-j shell, the highest angular momentum (denoted as I max ) 
state was found to have a sizable probability to be the g.s. fllP, [15|. This is explained 
by the observation that Afi max = 1 always: One can easily notice that the eigenvalue 
of / = Imax state is the lowest when Gj max = —1 and other parameters are switched 
off. Because Mj marL = 1, the predicted I max g.s. probabilities of fermions in a single-j 
shell are ~ = j^yo, x 100%) which are valid for all particle numbers (even or odd). 
It is predicted that the / = l ma x g-s. probabilities of fermions in a single-j shell 
decrease gradually with j and vanish at a large j limit; they will not saturate at a 
sizable value with j. 

Fig. 7a) shows the I m ax probabilities for different particle numbers in a single-j 
shell. The agreement between I max g.s. probabilities obtained by diagonalizing a 
TBRE hamiltonian and those predicted by using a simple jj x 100% is good. 

For (i-boson systems, the I ma x = 2n g.s. probabilities for all n are ~40-42%. In 
Sec. 2.1.2, the predicted P(/ ma:r )'s are A/j mox /5 = 40%, where Afi max = 2. 

For s<i-boson systems it was found in Ref. || that the I ma x g-s. probabilities 
are large, which can be actually explained in the same way. Among the two-body 
matrix elements, the interactions with C4 = — 1 and others being produce the 
lowest eigenvalue for the I max = 2n state. The predicted I = 2n g.s. probability 
is l/iV=l/6=16.7%, which is independent of the boson number. This is consistent 
with that obtained by diagonalizing a TBRE hamiltonian (~ 15%). Note that the 
term (s^S)(sd) gives degenerate lowest eigenvalues for many / states when e s d s d is 
set to be —1 and others are 0. Therefore, we use 6 (instead of 7) as the number 
of independent two-body matrix elements, N. The difference due to this minor 
modification is very small, though. 

For sdg-boson systems, the predicted l m ax = 4ng.s probabilities is 1/N ~ 3.2%, 
where iV = 32. The I max g.s. probabilities that we obtain by diagonalizing a TBRE 
hamiltonian are 3.3%, 4.2%, 3.3% for n =4, 5, 6, respectively. 

The above argument of P(/ max )'s can be generalized to more complicated cases, 
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such as fermions in a many-j shell, bosons with two or more different angular mo- 
menta (e.g., sdg bosons). Let us firstly take a two-j (j 1; j 2 ) shell. Similar to the 
above argument for fermions in a single- j shell, it is predicted that the two angular 
momenta I' max = I ma x(ji) and Imaxij^) have g.s. probabilities which are around or 
larger than 1/N x 100%. Here I max (j n ) refers to the maximum among all angular 
momenta of states constructed by j n configurations. In another word, one can pre- 
dict the lower limit of these I' max g.s. probabilities. Second, for a boson system, e.g., 
a s<ig-boson system, it is predicted that the / = I' max {d n ) = 2n g.s. probability is 
always larger than (or around) 1/32 x 100% = 3.2% (this probability obtained by 
diagonalizing a TBRE hamiltonian is ~12-13%). 

Fig. 7b) presents the I' max = I m axUi) , ImaxUZ) S- s - probabilities of fermions in 
a two-j shell, obtained by diagonalizing a TBRE hamiltonian. They are compared 
with the curve plotted using 1/N. It is noticed that the predicted lower limit of the 
{Imax)' g-s. probabilities works quite well. It is worthy to mention that other P(/)'s 
with / very near (I max )' are almost zero (smaller than 1%) in these examples. 

Now we study the eigenvalue of the I max state by calculating the aj max of fermions 
in a single- j shell. In doing so we present an argument of an observation that the 
a J T max is always lower than other eigenvalues of all other states while G j = — 1 
and others switched off. Equivalently, below we calculate aj^Z ^ setting Gj max =l 
and others zero, which gives the cyj^Z ^ ne largest among a/^'s of all J/3's. 

The calculation of otj™* is straightforward. By decoupling the two-body inter- 
action operators and using analytical formulas of Clebsch- Gordon coefficients, one 
can obtain all ai 's. It is noticed easily that there are only positive contributions 
in this state and there are always cancellations in all the other states when Gj max —1 
and others G/s switched off. The reason is as below: 

The wavefunction of the I max state is known as 

\ImaxM = I max ) = \jm 1 ,jm 2 , ■ ■ -jjmn) 

= ~ l)JU - 2), • • + 1 - n)). 

All / 7^ Imax states can be constructed by a successive orthogonalization with those 
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obtained by acting J_ operator on \I max M) state. It is easy to realize that both 
negative sign and positive sign appear in the wavefunctions of / 7^ Imax states. The 
coefficients in \I max M) may be chosen to be positive for all M. The aj max is given 
by a summation of squares of all possible couplings in the wavef unction. Therefore, 
there is no cancellation in calculating aj max . Cancellation appears if / 7^ I ma x- 

Below we list some results for n = 3 to 6 fermions in a single- j shell. 

1). For n = 3: 



a] 



2j 



a", 



2). For n = 4: 



a: 



Jmax 2 



= 3 + 



a 



-* max 



a' 1 



= 2 



2(4j-3)' 
3(2j - 2) 
2(4j-3)- 



2j(10j ~ 11) 
2(4j-3)(4j-5)' 

2(2j) 3 -11(2j) 2 + 9(2j) + 15 



(4j-3)(4j-5)(4j-7) 
5(2j-3)(2j-4) 
2(4j-5)(4j-7)' 

2j(512j 3 - 2848j 2 + 5116j - 2990) 



3 + 



8(8j-13)(4j-3)(4j-5)(4j-7) 



3). For n = 5: 



a, 



4 + 



4). For n = 6: 



a. 



(2j)(8-2j-17) 
2(4j-3)(4j-5) 
5(2j)(2j-1)(2j-2) 
8(4j-3)(4j-5)(4j-7)' 
35(2j-4)(2j-5)(2j-6) 
8(4j-7)(4j-9)(4j-11)' 



3 + 2j-3 + 148j 2 - 242j + 60 



4j 



4(4j-3)(4j-5) 



2j(2j - l)(396j 2 - 1482j + 1356) 
'8(4j-3)(4j-5)(4j-7)(4j-9)- 



(16) 



(17) 



(18) 



(19) 



21 



7 CORRELATION BETWEEN THE GROUND 
STATES OF FERMIONS IN A SINGLE-j SHELL 



It is interesting to see whether there is any correlation between different systems by 
the same two-body random interactions. For this purpose the cases of fermions in 
a single- j shell is interesting because, as shown in Sec. 2.2, despite of the simplicity 
these cases exhibit most regularities of many-body systems interacting via a TBRE. 

There are some evidences of correlation between ground states of systems with 
different particle number n but the same set of random interactions. In this section 
we firstly report a correlation of I g.s. probabilities for fermions in a single- j shell, 
then examine another correlation discussed in |B||, where this latter correlation was 



explained as reminiscence of the general seniority scheme |20 1 . 

The g.s. probability of 4 fermions in a single- j shell was found to fluctuate 



periodically with Sj = 3 |T(], [T5[]. We have checked n = 4 up to j =33/2, n = 6 up 
to j = 27/2 and n = 8 up to j = 23/2, which show that the g.s. probabilities 
change synchronously ||. The I = j g.s. probability of 5 fermions in a single- j 



shell exhibits a similar pattern |]16| . In this paper we notice that this synchronous 
fluctuation appear at an interval 5j = 3 for g.s. probability of systems with n = 4, 
6, and 8 fermions and I = j g.s. probability of odd n(=5,7) fermions in a single- j 
shell, but the origin of this correlation is not yet been available. 

In a previous work |J another correlation for fermions in a many-j shell was 
reported: the pairing phenomenon seems to be favored simply as a consequence of 
the two-body nature of the interaction. The "pairing" here means that there is a 
large matrix element between the S pair annihilation operator between the ground 
states of a n fermion system to an — 2, n — A ■ ■ ■ system. Below we examine this 
"pairing" correlation for fermions in a single-j shell, where the seniority quantum 
number v is well defined. 

First we see the simplest case: 4 and 6 fermions in the j = 11/2 shell. The 
g.s. probability for n = 4 and 6 is 41.2% and 66.4%, respectively. Among 1000 sets 
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of the TBRE hamiltonian, 364 sets give g.s. for both n = 4 and 6 simultaneously. 
This means that a TBRE hamiltonian in which I = is the ground state for n = 4 
has a extremely large probability (around 90%) produces 1 = ground state for 
n = 6. The overlaps between the g.s. of n = 6 and the g.s. of n = 4 coupled with 
an S pair for the same TBRE hamiltonian are in most cases around 0.8-0.9, while 
those between the g.s. of n = 4 and that of S pair acting on the g.s. of n = 6 
are in most cases around 0.9-1.0. This strongly supports the finding in ||: the 5* 
pair annihilation operator takes the g.s. of n fermions to that of n — 2 fermions. It 
is noted that the expectation values of seniority for g.s. is, more or less, randomly 
distributed from to 4, in these cases. Namely, there seems no bias for very low 
seniority in the above calculation by using a TBRE hamiltonian. 

By examining the g.s. of 4 and 6 fermions in the j = 11/2 shell, one notices 
that the seniority f's of these g.s. are quite close. For n = 6 there are strong 
seniority mixings between states with seniority v =0 and 4, but no mixing between 
states with v = 6 and (or 4). It is interesting to note that no TBRE hamiltonian 
of n = 6 for g.s. with seniority 6 produces g.s. of n = 4. It is noticed that the 
1 = state with v = 6 contributes ~ 24% to the total g.s. probability (66.4%) for 
n = 6. Roughly speaking, in this small single- j shell most of the random two-body 
interactions which produce g.s. with seniority v ~0-4 of n = 6 give g.s. with a 
similar seniority for n = 4, and in these cases the picture given in is appropriate. 

As for the j = 13/2 shell where there are strong mixings between states with 
v = 6 and v =0 (or 4), g.s. probability is 22.3% for n = 4 and 42.4% for n = 6. 
Among 1000 runs we obtain 13% sets of random interactions which give g.s. for 
both n = 4 and n = 6. Among those cases which simultaneously produce g.s. for 
n = 4 and 6, the overlaps between the g.s. of n = 6 and the g.s. of n = 4 coupled 
with an S pair for the same TBRE hamiltonian are in most cases around 0.6-0.9, 
while those between the g.s. of n = 4 and that of S pair acting on the g.s. of 
n = 6 are in most cases around 0.8-1.0, indicating a similar picture described in ||. 

Now we come to larger shells such as j = 15/2 to 23/2 with n =4, 6, and 8. 
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We concentrate on the j = 15/2 shell which is enough to demonstrate our points 
of view. The g.s. probability is 50.2%, 68.2% and 32.1%, for n = 4, 6 and 8, 
respectively. We have 31% (among the 1000 runs) of the TBRE parameters which 
produce g.s. for all n =4, 6 and 8, i.e., almost all those TBRE parameters which 
produce g.s. for n = 8 produce g.s. for n = 4 and 6. It is interesting to note 
that the expectation value of seniority v of all 31% g.s. of n = 4, 6 and 8 systems 
has a very small probability larger than 4, indicating a similar pattern given in [[|. 

On the other hand, it is noted that the the fluctuation of g.s. probabilities for 
n = 4, 6, and 8 fermions in the same single- j shell can be large. For example, the 
g.s. of 6 fermions in a j=15/2 shell is 68.2% while that of 8 fermions is 32.1%, 
which means that more than 50% of the g.s. for n = 6 are not related to the chain 
in which the g.s. of n — 2 fermions can be obtained by annihilating one S pair 
on that of n fermions. Fig. 8 presents a few cases of 4 or 6 fermions in a single- j 
shell, where no bias of a low seniority is observed, indicating that the contribution 
to the total g.s. probability beyond the a seniority chain described in || may be 
more important in realistic systems. For example, there are only 14% of the TBRE 
hamiltonian parameters which produce g.s. for both n = 4 and 6 fermions in a 
7=25/2, while the g.s. probability of n = 6 is 57.1%. 

Therefore, we conclude that the finding H of a chain of g.s. which is linked by 
S pair for fermions in a many-j shell is also observed frequently for even fermions 
in a single- j shell, and that this chain covers, however, only one part of the g.s. 
and the contribution beyond this chain may be more important. 

It is noted that no bias of low seniority in the g.s. is observed in our calculations 
by using a TBRE hamiltonian except that n = 8 fermions in a j = 15/2 shell, where 
no states with seniority 8 are observed in the 1 = ground states, i.e., most of 
g.s. have expectation value of v from to 4. 
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8 SUMMARY AND DISCUSSION 



To summarize, we have presented in this paper our understanding of regularities of 
many-body systems interacting via a two-body random ensemble. 

First, beginning with simple systems in which the relation between eigenvalues 
and two-body matrix elements is linear, we propose an integral to predict the I g.s. 
probabilities, P(J)'s. The properties of P(/)'s are understood by the shifts, or the 
largest and/or smallest eigenvalues with only one of the two-body matrix elements 
switched on and others switched off. This argument is further developed to predict 
the P(/)'s by a simple formula P(I) = Afi/N m , where A/} is the number of a certain 
angular momentum I gives the lowest eigenvalue among all the possible eigenvalues 
with one of the two-body matrix elements being —1 and others being zero, and N rn 
is taken as the number of two-body matrix elements except for very few cases. This 
method, as we show by a variety of very different systems, is applicable to both 
(even or odd number of) fermion systems (with both a single- j shell and a many-j 
shell) and boson systems. The agreement between the predicted I g.s. probabilities 
and those obtained by diagonalizing a TBRE hamiltonian is good. It is noted that 
this method predicts the g.s. probability, and on the same footing it addresses 
other I g.s. probabilities as well. Therefore, we provide in this paper a universal 
approach in studying the / g.s. probabilities. 

Next, we discuss the microscopic foundation of our simple approach by defining 
a local space defined by Gj = — 1 + {eGji(J' ^ J)}. We show that the angular 
momentum / which gives the lowest eigenvalue when Gj = — 1 and others are zero 
continue to be the lowest if e is small. One has the / g.s. probability around 70-90% 
even e is quite large (such as 0.5-1.0). This could be a naive explanation of the 
success of our approach, but a sound explanation is not yet available. 

A discovery of this work is that we are able to tell (by numerical experiments) 
which interactions, not only the monopole pairing interaction, are essential in fa- 
voring g.s. in both boson fermions systems. For instance, interactions with 
J = 0, 6, 8, 12, and 22 give the g.s. dominance for 4 fermions in a j = y shell. 
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This disproves a popular idea that the g.s. dominance essentially comes from the 
two-body nature of the interactions and might be independent of the form of the 
hamiltonian. 

The simple I ma x g-s. probabilities of fermions in a single-j shell are found, for the 
first time, to be determined by only the number of the two-body matrix elements 
and independent of particle number, and to follow a simply 1/N relation. This 
phenomenon is explained by the fact that Mi max = 1 for fermions in a single-j shell. 
A generalization of this regularity to fermions in a many-j shell and (sd and sdg) 
boson systems is shown to work well, too. Quite a few counter examples of the 
g.s. dominance in both fermions in a single-j shell or a multi-j shell and bosons are 
found and explained. 

One interesting note of the Tables V-IX is on quadruple pairing interaction for 
fermions in a single-j shell. It is found that I = n (j + 1 — n/2) give the lowest 
eigenvalue when G<i = — 1 and Gj^ — for all even (odd) number of fermions in 
a single-j shell that we have checked in this paper. It is not known whether this 
observation is always correct, and what the origin might be if it is. 

In this work we also studied the seniority distribution for fermions in a single-j 
shell where the seniority number is well defined. In a pioneering work || it was 
claimed that the (low) seniority chain is very important in the g.s. dominance, 
namely, the g.s. of n fermions is related to g.s. of n — 2, n — 4 and etc fermions via 
S pairs approximately. We confirm this phenomenon in this work. But we also note 
that the g.s. beyond this (approximate) seniority chain may be more important. 

It has not yet understood at a more microscopic level that why the above A/"o is 
large for even fermion systems, or even more specifically, why there is a staggering 
on the A/o for even fermions in a single-j shell. Further consideration of this issue is 
warranted. 

Finally, it is stressed that the P(/)'s discussed in this paper and the P(/)'s 
[ [Hfl (the probabilities of energies averaged over all the states for a fixed angular 
momentum I being the ground states) are different quantities. For even systems 
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the behavior of these two are accidentally similar. For odd- A systems, however, the 
■p(/)'s, which were explained in terms of the randomness of two-body cfp's, are very 
different from P(/)'s. This explicitly demonstrates that the / g.s. probabilities (and 
g.s. dominance) are not consequences of geometric chaoticity ||10|| . We show in 
this work that the g.s. dominance is actually related to two-body matrix elements 
which have specific features. 
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APPENDIX A Many-body systems interacting via a displaced TBRE 
While a TBRE is symmetric to zero, an interesting question is that what the 
results are if one uses random interactions with only positive, or negative sign, or 
random interactions which are not symmetric to zero. This issue is both interesting 
and important because interactions in realistic systems, such as nuclei, atoms etc, 
are not symmetric to zero. Below we present preliminary results by using a displaced 
TBRE. 

Let us firstly consider two arbitrary ensembles {G'j} and {Gj}, which are related 
by a shift C: 

G'j = Gj + C, J = 0,2,---2j-l, (20) 

where C is a constant. 

For fermions in a single- j shell, the results using the ensemble {G'j} are exactly 
the same as those obtained using {G'j} except a shift ^flc on the eigen-energy of 
the ground state. Therefore, a displacement of a TBRE is trivial in a single- j shell. 

For fermions in a two-j shell, however, the role played by the displacement of 
the TBRE is very complicated. Below we mention only two extreme situations: 

1. 4 fermions in a (2j 1} 2j 2 ) = (11, 3) shell, by using a TBRE, a TBRE+5, and a 
TBRE-5, the g.s. probabilities are 60.8%, 95.4%, and 100.0%, respectively. Here 
the g.s. is more pronounced if the TBRE is displaced to either negative or positive. 

2. 4 fermions in a (2j 1 ,2j 2 ) = (13,9) shell, by using a TBRE, a TBRE+5, and 
a TBRE-5, the g.s. probabilities are 44.8%, 2.7%, and 0.1%, respectively. Here 
the g.s. is greatly quenched down (close to zero) if the TBRE is displaced to 
either negative or positive. One observes a similar situation for 4 fermions in a 
(2j!,2j 2 ) = (7,5) shell. 

Concerning the effect of the shape of an ensemble, we consider the TBRE and 
an ensemble of uniformly changed random numbers between -1 and 1. The general 
features are quite similar. Slight differences appear if one multiply a factor to each 
Gj, such as the RQE of [2]. There seems no essential difference between them, 
however. 
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APPENDIX B Calculation of oi J Isg for 4 fermions in a single- j shell 



In this appendix, we calculate a™ for 4 fermions in a single- j shell by using 9j 



coefficients and coupled pair basis. 

The normalized pair basis of a 4-fermion system is defined as follows: 

\UL 2 : I) N = . 1 \A L ^ x A L ^Y |0> = . 1 [LiL 2 : J), (21) 



where the subscript N means that the state is normalized, and Nl 1 l 2 j is the overlap 
(LiL 2 : I\L\L 2 : I), which is given as follows: 

N LlL2l = l + 5 LlL2 -A{2L 1 + l) {2L 2 + 1)\ j j L 2 \ . (22) 

[ la I* I ) 

The matrix elements of two-body interaction within normalized pair basis are 
given in Eq. (A.26) of Ref. |22|: 

o , , 1 



<(L 1 L 2 :I\V2jTl[A J ^xA J \ \L 1 L 2 :I) n = ^— ]T {U LlL2 j R f (23) 

with 



NliL 2 I _R=even 



/ [00 


0)o = 


|00 


0)n, v = 0, 


\ |/o 


I)o = 


70 


I)n, v = 2 



~ \ 3 3 L1 \ 

U Li l 2 jr = 5 JLi Srl 2 + {-ySjL^RL, ~ ±L X L 2 JR I j j L 2 \. (24) 

{ J R I \ 

Here L x is a short hand notation of y/2L x + 1. 

The seniority v =0 and 2 states are easy to construct: 

(25) 

The subscript means the state is a seniority-conserved state. The seniority v — 4 
state is given by 

| ili2 : „„ = 1 | ili2 : _ « | /0 : 1),, (26) 
v 1 — « v 1 — or 

where Li 7^ 0, L 2 7^ 0, and a =n (10 : I\L\L 2 : Suppose there are more 

than one seniority 4 states, say, \L' X L' 2 : J)o and \L\L 2 : 7}o, one should or- 

thonormalize them to have an orthonormalized basis. We note here that in this 
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case the matrix elements q(L\L2 : I\y2J + 1 

o 



|LiL 2 : I) and (^i^2 : 



o are not an invariant, their summation is an in- 
variant. 

The q:q of seniority state is as follows: 



o(00 : 0|V2J + 1 U Jt x A J ] |00 : 0) 



8 



o 2 -?- 1 7 = 

~ 2 i+l 



4j -2_ 1; 



o 



(27) 



The 9j coefficients are difficult to be further simplified unless one of them are 0. 
The ai max is derived in another way in Sect. III. 4. 
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APPENDIX C A few interesting examples 



The g.s. dominance was reported to be robust in even fermion systems and 
s<i-boson systems. Similarly, for fermions in a single- j shell it was believed that the 



I — j g.s. probability is the largest |10| . However, as shown in simple cases such as 
c?-boson systems, these dominances are not always true. Here we list a few counter 
examples that we noticed in our work. 

1. 4 fermions in a j — | shell, which was explained in terms of shifts in fllS] ; and 
4 fermions in a j = y shell, which is understood in this paper by counting Aq. In 
these two examples, P(2)'s are larger than -P(O). 

2. A two-j shell with (2ji,2j2) = (7,5) and n = 4, where the g.s. probability 
is 21.5% and that of 2 g.s. is 34.7%; 

3. two-j shells with (2j 1? 2j 2 ) = (H,7), (13,9) and n = 5, where ji g.s. and j 2 
g.s. probabilities are small (< 20%) while | g.s. probability is large (> 30%), which 
means that an angular momentum which is not ji or j 2 may be favored to be the 
ground state for odd number of fermion systems; 

4. d-bosons systems with n ^ Qn, while the g.s. probability is less than 40% 
(2 % or ~ 36%) while the I max g.s. probability is ~ 42%. Especially, the cases 
with n = Gk ± 1 (k, is a natural number) are extreme counter examples of the 
g.s. dominance: the g.s. probabilities are very close to zero periodically. All 
regularities of -P(/)'s of <i-boson systems are well understood in this paper. 
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TABLE I. The probability of each state to be the ground state and distribution 
width of each eigen-energy in the case of j = 5/2 shell with 3 fermions. All the 
states are labeled uniquely by their angular momenta /. The probabilities of the 
row "TBRE" are obtained by 1000 runs of a TBRE hamiltonian, and those of 
"predl." are obtained by integrals such as Eq. (0) for + state of n = 4,j = | case. 
The row "pred2." are obtained by the new approach proposed in this paper. The 



distribution width, gj, of each eigen-energy, is listed in the last row. 



/ 


3/2 


5/2 


9/2 






TBRE 


40.1% 


23.7% 


36.2% 






predl. 


41.82% 


22.77% 


36.37% 






pred2. 


40% 


20% 


40% 






9i 


2.31 


1.84 


2.44 






TABLE II. Same 


as Table I except that j = 


7/2 shell. 




I 


3/2 


5/2 


7/2 9/2 


11/2 


15/2 


TBRE 


30.4% 


29.7% 


9.6% 4.4% 


0.0% 


25.9% 


predl. 


31.14% 


27.26% 


9.70% 3.69% 


0.00% 


28.13% 


pred2. 


28.6% 


28.6% 


14.3% 





28.6% 


9i 


2.44 


2.09 


1.57 2.11 


1.88 


2.42 



TABLE III. Same at Table II except that n — 4. All the eigenstates are uniquely 



labeled by I(v). 



I(v) 


0(0) 


2(2) 


2(4) 


4(2) 


4(4) 


5(4) 


6(2) 


8(4) 


TBRE 


19.9% 


1.2% 


31.7% 


0.0% 


25.0% 


0.0% 


0.0% 


22.2% 


predl. 


18.19% 


0.89% 


33.25% 


0.00% 


22.96% 


0.00% 


0.02% 


24.15% 


pred2. 


14.3% 


0% 


28.6% 


0% 


28.6% 





0% 


28.6% 


9i(v) 


3.14 


3.25 


4.12 


3.45 


3.68 


3.62 


3.64 


4.22 
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TABLE IV. The angular momenta which give the largest (smallest) eigenvalues 
when c\ — — 1 and other parameters are for d boson systems. 



n 


c (mm J 


C2 (mmj 


c 2 (maxj 


C4(mm 


6k 








Imax 


Iraax 


6k+1 


2 


2 


Imax 


Imax 


6k+2 





2 


Imax 


Imax 


6k+3 


2 





Imax 


Iraax 


6k+A 





2 


Imax 


Imax 


Qk+5 


2 


2 


Imax 


Imax 



cA max) 




2 
2 

2 

9 



TABLE V. The angular momenta which give the lowest eigenvalues when Gj = 
— 1 and other parameters are for a 3-nucleon system in a single- j shell. Here we 
use a unit of 21, and M refers to I max - 



2j 


Go 


G 2 




G 4 


G 6 




G 8 


Gio 


Gl2 


G14 


Gl6 


Gig 


G20 


G22 


G24 


G26 


G28 


G30 


5 


2j 


2j - 


2 


M 






























7 


2j 


2j - 


2 


3 


M 




























9 


2j 


2j - 


2 


2j - 6 


9 




M 
























11 


2j 


2j - 


2 


2j - 6 


3 




15 


M 






















13 


2j 


23 - 


2 


2j - 6 


3 




9 


21 


M 




















15 


2j 


23 - 


2 


2j - 2 


5 




3 


15 


27 


M 


















17 


2j 


2j - 


2 


2j - 2 


2j - 


6 


3 


9 


17 


33 


M 
















19 


2j 


2j - 


2 


2j - 2 


2i - 


6 


5 


3 


15 


23 


M-6 


M 














21 


2j 


2j - 


2 


2j - 2 


2j - 


6 


11 


3 


9 


17 


29 


M-6 


M 












23 


2j 


2j - 


2 


2j - 2 


2j - 


6 


13 


5 


3 


11 


23 


35 


M-6 


M 










25 


2j 


23 - 


2 


2j - 2 


2i - 


6 


15 


11 


3 


9 


17 


29 


41 


M-6 


M 








27 


2j 


2j - 


2 


2j - 2 


2j - 


6 


29 


13 


5 


3 


11 


23 


35 


47 


M-6 


M 






29 


2j 


23 - 


2 


2j - 2 


2j - 


2 


2j - 6 


15 


11 


3 


9 


17 


29 


41 


53 


M-6 


M 




31 


2j 


2j - 


2 


2j - 2 


2j - 


2 


2j - 6 


21 


13 


5 


3 


11 


23 


35 


47 


59 


M-6 


M 



35 



TABLE VI. The angular momenta which give the lowest eigenvalues when Gj 
■1 and other parameters are for 4 fermions in a single- j shell. 



2j 


Go 


G 2 


G 4 


G 6 


G s 


Gio 


G12 


Gl4 


Gi6 


Gi8 


G20 


G22 


G24 


G 2 6 


G 2 s 


G30 


7 





4 


2 


8 


























9 





4 








12 
























11 





4 





4 


8 


16 






















13 





4 





2 


2 


12 


20 




















15 





4 





2 








16 


24 


















17 





4 


6 





4 


2 





20 


28 
















19 





4 


8 





2 


8 


2 


16 


24 


32 














21 





4 


8 





2 











20 


28 


36 












23 





4 


8 





2 





10 


2 





24 


32 


40 










25 





4 


8 





2 


4 


8 


10 


6 





28 


36 


44 








27 





4 


8 





2 


4 


2 








4 


20 


32 


40 


48 






29 





4 


8 








2 


6 


8 


12 


8 





24 


36 


44 


52 




31 





4 


8 








2 





8 


14 


16 


6 





32 


40 


48 


56 



TABLE VII. Same as TABLE V except that n = 5. 



2j 


Go 


G 2 




G4 




Gs 


C10 


C12 


G14 


Ci6 


Ci8 G20 


9 


2j 


2j- 


4 


9 


3 


Imax 












11 


2j 


2j- 


4 


5 


2j 


5 


Imax 










13 


2j 


2.7 - 


4 


5 


2j 


7 


5 


Imax 








15 


2j 


2 j - 


4 


7 


9 


2j 


7 


31 


Imax 






17 


2j 


2j- 


4 


9 


5 


2j 


11 


9 


41 


I-max 




19 


2j 


2j- 


4 


11 


5 


13 


2j 


13 


5 


51 


Imax 


21 


2j 


2j- 


4 


19 


5 


7 


2j 


15 


2j 


2j 


6 1 I max 
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TABLE VIII. Same as TABLE V except that n = 6. 

2j Gq G2 G4 Gq Gs Gio G12 G14 G\o Gi$ G20 



11 





6 


4 








Imax 








13 





6 


4 





4 


4 


Imax 






15 





6 





6 











Imax 




17 





6 











4 





16 


Imax 


19 





6 





6 











4 


22 T 

1 max 


21 





6 


10 


6 
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TABLE IX. Same as TABLE V except that n = 7. 



2j 


Go 


G 2 




G 4 


Go 


G 8 


C10 


C12 


G14 


Gi6 Gi8 


13 


2j 


2j- 


6 


11 


1 


11 


5 


Imax 






15 


2j 


2j- 


6 


13 


3 


3 


2j 


9 


Imax 




17 


2j 


2j- 


6 


7 


2j 


2j 


13 


9 


23 


Imax 



TABLE X. Same as TABLE V except that s<i-boson systems. The matrix element 
corresponding to e s d s d are omitted because it always presents degenerate levels for 
ma ny I states. 



n 


&SSSS 


&sddd 


&ssdd 


c 


c 2 


c 4 


6 

















Ifnax 


7 











2 


2 


Imax 


8 














2 


Imax 


9 











2 





Imax 


10 














2 


Imax 


11 











2 


2 


Imax 


12 

















Imax 


13 











2 


2 


Imax 


14 














2 


Imax 


15 











2 





Imax 


16 














2 


Imax 
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Figure captions: 

FIG. 1 The I g.s. probabilities of d bosons. The boson number n runs from 4 
to 44. Only states with / =0, 2, and I max — % n are possible to be the ground. The 
g.s., 2 g.s. and I max = 2n g.s. probabilities are very near to 0, 20%, 40% or 60%. The 
P(0) ~ in rf-boson systems with = 6k ± 1. The predicted P(/)'s (open squares) 
are well consistent with those (solid squares) obtained by diagonalizing a TBRE 
hamiltonian. All regularities are explained by the reduction rule of U(5) — > 0(3). 

FIG. 2 The P(/)'s of fermions in a single-j shell. The solid squares are P(/)'s 
obtained by diagonalizing a TBRE hamiltonian and the open squares are P(/)'s 
predicted by the approach proposed in this paper. a) j — | with 5 fermions; b) 
j — | with 4 fermions. Good agreements are obtained for both even and odd-A 
cases. 

FIG. 3 The P(0)'s of fermions in a single-j shell. Solid squares are obtained 
from 1000 runs of a TBRE hamiltonian. The open squares are predicted P(0)'s in 
this paper. a) n = 4; b) n — 6. solid triangles in a) are obtained from an 
empirical formula, Eq. (15). 

FIG. 4 The P(0), P(2) and P(I max ) of srf-boson systems. Solid symbols are 
P(/)'s obtained from 1000 runs of a TBRE hamiltonian. Open symbols are P(/)'s 
predicted in this paper. Only / = 0, 2, I max g-S- probabilities are included. All 
other P(/)'s obtained by diagonalizing a TBRE hamiltonian are close to zero, and 
the predicted P(/)'s are zero. 

FIG. 5 Fermions in a two-j shell with (ji, j 2 ) = (|, f). n =4, 5, 6, 7 in a), 
b), c) and d), respectively. Solid squares are obtained by 1000 runs of a TBRE 
hamiltonian and open squares are predicted by the method proposed in this paper. 

FIG. 6 4 fermions in a j = y shell. a) The I max g.s. probabilities with 
Gi 6 =±l (J max =^Q) and all other Gj being a TBRE multiplied by e; b) The g.s. 
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probabilities with G =±l and all other G j being a TBRE multiplied by e. Refer to 
text for details. 

FIG. 7 The I max g.s. probabilities, a) The I max g.s. probabilities of fermions 
in a single- j shell, which are obtained by diagonalizing a TBRE hamiltonian (1000 
runs). Solid squares are plotted by x 100%, N — j + |. It is noted that the I max 
g.s. probabilities follow a x 100% relation well, and are independent of particle 
number, n. b) The (I ma x)' = Imax{j\) and Imaxij^) g- s - probabilities of fermions 
in two-j shells, obtained by 1000 runs of a TBRE hamiltonian. The lower limit of 
the (Imax)' g-s- probabilities are predicted x 100% (solid squares). The I' max g.s. 
probabilities are reasonably consistent with the predictions. 

FIG. 8 The seniority distribution in the angular momentum 1 = ground 
states. No bias of low seniority is observed in these 4 and 6 fermions in a single-j 
shell, which indicates that the contribution to the total g.s. beyond a low seniority 
chain may be more important. 
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